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Abstract 

The canonical correlation or Kubo-Mori scalar product on the state space of a finite quantum 
system is a natural generalization of the classical Fisher metric. This metric is induced by the 
von Neumann entropy or the relative entropy of the quantum mechanical states. An important 
conjecture of Petz that the scalar curvature of the state space with Kubo-Mori scalar product 
as Riemannian metric is monotone with respect to the majorisation relation of states: the scalar 
curvature is increases if one goes to more mixed states. We give an appropriate grouping for 
the summands in the expression for the scalar curvature. The conjecture will follows from the 
monotonicity of the summands. We prove the monotonicity for some of these summands and we 
give numerical evidences that the remaining terms are monotone too. Note that the real density 
matrices form a submanifold of the complex density matrices. We prove that if Petz's conjecture 
true for complex density matrices then it is true for real density matrices too. 

1 Introduction 

The state space of a finite quantum system can be endowed with a differentiable structure |l(Jj . The 
canonical correlation defines a Riemannian structure on it. There is strong connection between the 
scalar curvature of this manifold at a given state and statistical distinguishability and uncertainty of 
the state Roughly speaking the scalar curvature measures the average statistical uncertainty. 

This idea comes from a series expansion of the volume of the geodesic ball. If Dq is a given point in 
the state space and B r (Do) is the geodesic ball with center Dq and radius r then the volume of this 
ball is given by 

V r\ OJJ r (| + 1 ) ^ 6(71 + 2) K ') y ' 

where Scal(Do) is the scalar curvature at the point Dq and n is the dimension of the manifold. There 
was given an explicit formula for the scalar curvature of these manifolds, for example in [3J [7J |!|] . 
Physically it is reasonable to expect that the most mixed states are less distinguishable from the 
neighboring ones than the less mixed states, for details see |l()l 112) . It means mathematically that 
the scalar curvature of physically relevant Riemann structures should have monotonicity property in 
the sense that if D\ is more mixed than D2 then Scal(Da) should be less then Scal(Di). This was 
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conjectured first by Petz ^Uj and it was verified for 2x2 matrices in JU]. There were some numerical 
simulations showing that the conjecture is true but no mathematical proof for that. The aim of the 
paper is to give an appropriate grouping for the summands in the expression for the scalar curvature 
and prove that some of these summands are monotone with respect to the majorisation. 

The paper is organized as follows. In Section 2 definition, properties and equivalent forms of 
majorisation are given and the conjecture is formulated. In Section 3 the Kubo-Mori Riemannian 
inner product is defined and the curvature formula is given for real and complex density matrix 
spaces. In Section 4 we rewrite the scalar curvature formula as the sum of five terms. Four of them 
correspond with the conjecture as elementary, but brutal computations show. Further remarks and 
a lemma about the conjecture is presented as well. In Section 5 we prove that if the scalar curvature 
is monotone with respect to the majorisation on the space of complex density matrices then it is 
monotone on the space of real ones. 



2 Majorisation 

The state space of a finite quantum system is the set of positive semidefinite n x n matrices of trace 
1. Such matrices are often called density matrices. Let denote the space of invertible density 
matrices. The state A is called majorised by the state B, denoted by A -< B if the following hold for 
their decreasingly ordered set of eigenvalues (a\, . . . , a n ) and (b%, . . . , b n ) 

k k 

X><X> (2) 

1=1 i=i 

for all 1 < k < n. For any (ai, . . . , a n ) eigenvalues of a density matrix 



1 1 
, . . . , 

n n 



-< (ai,...,o„) -< (1,0,. ..,0). 



The majorisation occurs naturally in various contexts. For example, in statistical quantum mechanics 
the relation A -< B is interpreted to mean that the state A describes a more mixed or a " more chaotic" 
state than B £Q. A linear map T on R™ a T-transform if there exists < t < 1 and indices k,l such 
that 

T(x x , ...,x n ) = (xi, . . .,x k -i,tx k + (1 - t)xi,x k +i, ■ ■ - ,xi-i, (1 - t)x k + tx t ,xi + i, . . . ,x n ). 

Note that if the density matrix A has eigenvalues (a± , . . . , a n ) and the density matrix B has eigenvalues 
T(a±, . . . , a n ) then B < A. 

From every self-adjoint H operator and positive parameter (3 one can set a state 

Rntf) = (3) 

which is called the Gibbs state at the inverse temperature (3 for the Hamiltonian H. One can check 
that if the difference Hi — H2 is a multiple of the identity then Rh 1 (f3) = Rh 2 (/3). It is known that if 
Pi < P2 then Rh{Pi) ~< Rh{Pi) 8-,. This last result means that the states are more mixed at higher 
temperature. 

Theorem 2.1. Assume that we have two invertible states A and B with decreasingly ordered set of 
eigenvalues (oj_, . . . , a n ) and (pi, . . . , b n ). The following are equivalent: 

1. The state A is more mixed than B. 



2 



2. One can find a sequence (C z )z=l,...,d between them such that for all z = 1, . . . , d: C z € Q^ , 

A = C 1 <C 2 <---<C d = B 
holds and the set of eigenvalues of C z and C z -i is the same except two elements. 

3. The set [ax, . . . , a n ) is obtained from (&i, . . . , b n ) by a finite number of T-transforms. 
4- There is a sequence (G z ) z= i d between them such that for all z = 1, . . . , d: G z £ Q+, 

A = G x -< G 2 -< <G d = B 

holds and for all i = 1, . . . , d — 1 there exists a selfadjoint operator Hi and positive parameters 
Pi,uth,i such that d = RH t (Jh,i) and G i+ i = Rh,{P2,i)- 

Proof. That (1) and (3) are equivalent can be found for example in ^2|- The statement (2) is just the 
reformulation of (3). The implication (3) — > (2) is trivial and easy to check that (2) — > (3) implication 
holds too. □ 

After this introduction we can formulate correctly the Petz's conjecture |10j . 

Conjecture 2.1. If D\ -< D 2 then Scal{D\) > Scal(D 2 ) where Scal(D) denotes the scalar curvature 
of the state space induced by the canonical correlation inner product as Riemannian metric at the 
point D. 

3 Scalar curvature formula 

Let be the space of all complex self-adjoint positive definite n x n matrices and let M. n be 
the real vector space of all self-adjoint n x n matrices. The space JA^ can be endowed with a 
differentiable structure and the tangent space T^M.^ at D £ M.„ can be identified with M n [E] 
One can consider the quantum mechanical state space as a Riemannian submanifold of Ai^ of 
codimension 1. 

There is a very important functional on the space M.„, namely the von Neumann entropy 

S(D) = -TrDlog(D). (4) 
Since this functional is strictly concave [HIE]) the second derivative of the entropy 



ddS:M+^>UN(M n ><M n> R) D^((X,Y)^-J^ Tr((D + t)~ 1 X{D + t)~ 1 Y)^J 



(5) 



is negative definite. The Riemannian metric which arises as the negative of the second derivative of 
the entropy 

poo 

(6) 



POO 

G D (X,Y)= Tr((I3 + t)- 1 X(£) + t)- 1 y)df 
Jq 



is called Kubo-Mori metric. This scalar product is an important ingredient of linear response theory 
and often called the canonical correlation of X and Y. 

Let us introduce this metric in a different way. Recall that Umegaki's relative entropy 

S(Dx > D 2 )=TrD 1 (log D 1 - log D 2 ) (7) 
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of density matrices measures the information between the corresponding states jS] . The partial deriva- 
tives of the relative entropy is the Kubo-Mori inner product: 

d 2 

S(D + tX,D + sY) =G D {X,Y). (8) 



dtds 



\t=s= 



Let D € M.+ and choose a basis of W 1 such that D = J2k=i ^kEkk is diagonal, where (-Ey&)i,fc=i,...,n 
are the usual system of matrix units. Define the self-adjoint matrices 

Fm = Eui + Eik Hki = i Eki — i Eik , 

then the {F k i)i<k<i< n , {Hkl)x<k<l<n and (F k k)i<k<n vector system is a basis in M n . Let us define 
for positive numbers x, y, z the functions 

f°° 1 f°° 1 
m(x.y) := / ^ — -dt m(x,y,z):= / — — -dt. (9) 

V ,y> J (x + t)(y + t) 1 ' y ' j Jo + + + V ' 

The scalar product of the basis vectors in the tangent space at the point D: 

r 1 ^ ■ ^ • ^ 1 ^ 1 ^ 1 ^ j Go(E[ij, Hki) = GniFij, Fki) — SikSji2m(Xi, A,) 

tor 1 < s < j < n, 1 < k < I < n : < ) J / VJ J ■" 

y UD{tlij, rkl) — U 

for 1 < i < j < n, l<k<n: G D (H tJ ,F kk ) = G(F tJ ,F kk ) = (10) 

for 1 < i < n, l<k<n: G D (F tl ,F k k) = 6 tk 4m(\ l ,\ t ) 

where 8ik = 1 if i = k else dik = 0. We will use the following properties of the m functions: assume 
that x,y,z and fi are different positive numbers, then 



logx-logy m(x,z) - m(y,z) 1 1 

m(x,y) — , m(x,y,z) = , m(x,x) = — m(x,x,x) = — 3 

x — y x — y x 2x z 

m(x,y) = -m((ix,[iy), m(x,x,y) = - — -, m,(x,y, z) = ^m((ix, (iy, fiz). 

\x x-y ^ 

The following identity will be used several times 

1 f m(x 1 x,y) m{x,y,y) 



Let us define two functions 



m(x,y) \ m(x,x) m{y,y) 



(12) 



/ 2 ) = 1 m{x,y,z) 2 m(y,y,x)m(y,y,z) 

2 m(x,y)m(y 7 z)m(z 7 x) m(y , x)m(y , y)m(y , z) ' 



v(x,y) 



1 m(x,x,y) 2 m(x,x,y)m{y,y,x) 



2 m(x 1 y) 2 m(x,x) m{x,y) 2 m{y,y) 
From the definition and the identities (|ll(l one can derive scaling properties 

</j(/ia;, /iy, //z) = -99(2;, y, z) v{fxx, fxy) = -v(x, y). (14) 



The space Q+ is a one codimensional submanifold of M^. The tangent space of Q+ is the set of 
self-adjoint traceless matrices. The scalar curvature of this submanifold can be computed using the 
Gauss equation 015]. 
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The scalar curvature the space of complex density matrices at a given matrix D with eigenvalues 
Ai, . . . , A n is 

n 

Scal(£>) = V(*j,*k,\i), (15) 

j,k,i=i 
\{j,k,i}\>i 

where \{j, k,l}\ > 1 means that the indices i, j and k are not equal [3], and that in the space of real 
density matrices is 

Scal R L> = - V(^,A fc ,A,) + 2 U ( A *' A 0- (16) 



IO'.*.i}|>l 



j,k,l—l k,l-l 



[7]- Note that the scalar curvature at a given point depends only on the eigenvalues of the density 
matrix. 

4 About monotonicity of the scalar curvature 

Using Theorem 12 . II the monotonicity conjecture will follows from the inequality 
Scal(A) = Scal(Ci) < Scal(C 2 ) < < Scal(C* d ) = Scal(S). 

Let A and B be two states, such that A -< B and the decreasingly ordered set of eigenvalues of 
A and B be the same except two elements. In this case their eigenvalues can be written in the form 
(Ai, . . . , a, . . . , b, . . . , A n ) and (Ai, . . . , a — x, . . . , b + x, . . . , A n ), where x < 

If one computes the scalar curvature using equation (|15f) the summands can be grouped: 

1. Summands where only a and b appear 

a(a, b) := 2<p(a, a, b) + 2(p(b, b, a) + <p(a, b, a) + fib, a, b). (17) 

2. Summands where a,b and another eigenvalue appears 

/?i,k(a, b) := 2<p(a, a, X k ) + 2cp(b, b, X k ) + <ys(a, X k , a) + <p{b, X k ,b) (18) 
/?2,fc(a, b) := 2cp(a, b, X k ) + 2cp(b, a, Xk) + <p{a, X k ,b) + ip{b, X k , a). (19) 



(20) 



3. Summands where a or b appear only once 

7*2(0,6) :=ip(a, X k , A;) +(p(b,X k ,Xi) +ip(X k ,a,Xi) 

+ <p(X k ,b, Xi) + ip(X k , A ; , a) + ip(X k ,X h b). 

4. Summands without a and b 

Sj k i := (f(Xj,X k ,Xi). 
After this grouping the curvature formula is 

n n n 

Scal(D) = a(a, b) + b) + /3 2>fe (a, b)) + £ j k i(a,b)+ £ 5 jM . (21) 

k ~ 1 k,l-l j,kd=l 

\{j,k,l}\>l 
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One possible way to prove the monotonicity conjecture is to show that every summand in the 
previous formula is monotone with respect to the majorisation. 

Dittmann used the symmetrization of function <p and the previous grouping of the scalar curva- 
ture in 0. Numerical tests were confirmed by Dittmann about Petz's monotonicity conjecture and 
suggestive 3D-plots were given in |3] . 



Theorem 4.1. Part a of the scalar curvature is monotone, that is for given a > b positive numbers 
the function a(a — x, b + x) is strictly increasing in the variable x G [0, ^rp] ■ 

Proof. From the identities l|ll[) . (|13[l it follows, that 

-1 1 + b/a 
2(a + b) m(l,b/a) 2 



<p(a,b,a) + ip(b,a,b) = — — — rr ^^vL Ml, l,b/a) 2 + (b/a)m(b/a,b/a, l) 2 ) 



and 



„ / „ /, , x ~ 1 / , / ,oTn(l,l,b/a)m(b/a,b/a,l) 
2ip(a,a,b) + 2(p(b,b,a) = -(1 + b/a) 2 v ' 



a + b n(l,b/a) 2 
From these formulas and the special values of the functions m in (jl 1|> we obtain 

-1 / l(l + c(z)) 2 (1 + c{x))(l + c{x) 2 ) 1 (l + c(a;)) 2 



a(a — x, b + x) 



a + b \ 2{l-c{x)) 2 c{x){c{x) - l)logc(aj) 2 c (x) log 2 c(x) ) ' 
where 

, \ b + x 

c{x) = . 

a — x 

Since the function c(x) increasing, to prove the theorem it is enough to show that the function 

n l (l±£f | (i + c)(i + c 2 ) i (i + c) 2 

[ ' 2(l-c) 2 + c(c- l)logc 2 C l g 2 c 
is decreasing on the interval I=[0, 1]. To show that the function 

(c+1) 2 (c+ l)(3c 2 - 2c + 3) c 4 - 2c 3 - 2c 2 - 2c+ 1 2c + 2 



r'(c) 



c 2 log 3 c c 2 (c- 1) log 2 c c 2 (c- l) 2 log c (c-1) 3 

is negative or equivalently the function 

n(c) = (c -l) 3 log 3 c-t'(c) 

is negative on the interval /, enough to prove that 

(1) t[(c) > on / and lim c ^i ti(c) = 0. 

The limit can be easily checked and the first part will follow from the statement: 

(2) r 2 (c) > on I and lim^i r{(c) = 0, 
where 

T2(C) = (c-l)( C 2 - C +l)(c 2 + C +l) ' T i (c) - 
Substituting r{ into the previous formula: 

, 2 (c-l)(c+l)(c 2 -12c+l) , (c+l) 2 (c-l) 2 

r 2 c = 61og 2 c + , 9 A TT 1 ■ lo 8' c + ^TT^ TW? 9 7v3 

w (c 2 + c+ l)(c 2 - c+ 1) 2(c 2 + c+ l) 2 (c 2 - c+ l) 2 
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The limit in the statement (2) again can be checked and the positivity of r 2 (c) will follow from the 
next statement 

(3) r 3 (c) < on / and lim^i t 2 (c) = 0, 
where 

r 3 (c) = c(c 2 + c + l) 2 (c 2 - c + l) 2 • r 2 (c). 

Computation shows that 

r 3 (c) =2(6c 8 + 6c 7 + 13c 6 - 24c 5 + 22c 4 - 24c 3 + 13c 2 + 6c + 6) logc 
+ (c 8 - 12c 7 + 4c 6 - 4c 2 + 12c - 1). 

The limit in the statement (3) again can be checked and the negativity of r 3 (c) will follow from the 
statement 

(4) r 3 (c) > on I and lim c ^i r 3 (c) = 0. 

Easy to check that the limit condition fulfils in the previous statement, and the inequality is the 
following: 



4(24c 7 + 21c 6 + 39c 5 - 60c 4 + 44c 3 - 36c 2 + 13c + 3) logc 
(r — I) 

+ 2^ ^(10c 7 - 26c 6 - c 5 - 25c 4 - 3c 3 - 27c 2 - 18c - 6) < 0. 



c 



One can show that the coefficient of logc is strictly positive on /: The function /i(c) = c(100c 2 — 
144c + 52) is positive on /. The function /2(c) = 7c 6 + 12c 5 — 14c 4 + 2 has two stationary points on /: 
a local maximum at the origin and a local minimum at c = ^ 15 +^^ . and /2(co) > 0. The coefficient 
of log c is 

96c 7 + 6/ 2 (c) + h(c) + 36c 3 (l - c) 

which is positive on /. 

Taking into account that the coefficients of log c is strictly positive on I one can rearrange the 
terms in the previous inequality: 

2(c-l)(10c 7 -26c 6 -c 5 -25c 4 -3c 3 -27c 2 - 18c - 6) 
q\c) - og(c) + c( - 96c 7 + 84c 6 + 156c 5 _ 240c 4 + 176c 3 - 144c 2 + 52c + 12) ~ ^ ' 

We will show this inequality in two steps. 

1. < c < 2 : In this case let decrease the function q(c) and show that the inequality 

-•(c) = log(c) + 3(c-l)(2c 7 -l) > 

y w 6 w c(24c 7 + 24c 6 + 44c 5 - 60c 4 + 44c 3 - 36c 2 + 13c + 13) ~ 

holds. Since q* (|) > 0, it is enough to show that dq Jf^ < if < c < |. This follows from the 
inequality 

dq*(c) _ -1 



dc c 2 (24c 7 + 24c 6 + 44c 5 - 60c 4 + 44c 3 - 36c 2 + 13c + 3) 2 

• (c 12 (-576c 3 - 1440c 2 - 3216c + 2112) + c 10 (-2944c + 1472) 

+ c 8 (5296c 2 - 7700c + 2640) + c 6 (4800c 2 - 7292c + 2508) 

+ c 3 (1800c 3 - 568c 2 - 624c + 241) + (550c 3 - 441c 2 + 69c + 9)) , 

where dq is the sum of six negative functions. 
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2. | < c < 1 : Since q(l) = enough to show that q'(c) < if | < c < 1. Let us define the 
following function 

k(c) = -c 2 (24c 7 + 24c 6 + 44c 5 - 60c 4 + 44c 3 - 36c 2 + 13c + 3) 2 • q'(c). 

The aim is to show that k(c) is positive if | < c < 1. Since 

Jl) .,(1) .^(^ «<»>m>o 



2) 1 \2J ' \2J ' 7 V 2 , 
enough to show that «( 10 )(c) >0if 1 j<c<l. This comes from the equality 

K ( 10 ) ( c ) = 39916800(104832c 4 + 110292c 3 + 50688c 2 + 3963c + 1015). 



□ 



Theorem 4.2. Part ^x,h of the scalar curvature is monotone, that is for given 1 > a > b positive 
numbers and eigenvalue 1 > Xk > the function /?i,fc(a — x, b + x) is strictly increasing in the variable 



Proof. From the identities (| 1 1 1) and (|13|l it follows, that 



2<p(a,a, X k ) + <p(a, A fe ,a) = ■ r(a/A fe ), 



where 

2c - 3 

r(c) 



2c log 2 c c(l-c)logc 2(1 -c) 2 ' 
From these equations one arrives at the expression 

n , , \ 1 f fa — x\ [b + x 
Pi,fc(a- x,b + x) = — ■ t 



Xk \ \ Xk J \ X 

This means that Pi.k(a — x, b + x) is increasing if r(a — x) + t(6 + x) is increasing on [0, 3 ^], where 
a = a/Xk, b = b/Xk and x — x/Xk- To prove that 

-r'(a - x) + r'(b + x) > 

holds, which means that r'(x) decreasing, we will show that 

r"(x) < 0. 

We will do this in two steps: we show that the functions 

Tl(C) = icbg^" 2c(c-l)logc' T2 (C) = 2ck^Tc ~ c(c - 1) log c + (T^) 2 (23) 

are concave, and we note that ri(c) + (1/2)t2(c) = t(c). 
In the concavity proof of t\ (c) we work with the function 

t* (c) = -2c 3 (c -l) 3 log 4 c-r{'(c) 

that is 

r*(c) = (6x 2 - 6x + 2)log 3 x- (x- l)(3x- 2)(x- 3) log 2 x 
+ (x - l) 2 (x 2 - lOx + 11) logx + 3(x - l) 3 (x - 3). 



To prove the concavity of t\ (c) enough to show that t* (c) is negative if c G / and positive if c > 1 . 
Since lim c ^i r*(c) = enough to show that r*W(c) > 0. 

Since lim c _>i r*'^(c) = enough to show that t*^ 2 \c) is negative if c S I and positive if c > 1. 
Since lim c ^i t*^(c) — enough to show that r*( 3 )(c) > 0, or equivalently 

p(c) =c 3 ■ r* (3) (c) = (-18c 3 + 36c 2 - 18c + 12) log 2 c 

+ (24c 4 - 138c 3 + 164c 2 + 34c - 12) log c + (98c 4 - 276c 3 - 184c 2 - 4c - 2) > 0. 

Using the previous method the positivity of p(c) comes from the positivity of the function p^(c) 
and the limits 

lim p (0) (c) = lim p (1) (c) = lim p (2) (c) = lim p (3) (c) = 0. 

c — >1 c — >1 c — >1 c — >1 

The limits can be checked. The positivity of 

p [i) (c) = -j (72(8c 4 - 3c 3 - 2c 2 + c - 2) log c + 8(444c 4 - 153c 3 - 50c 2 + 4c + 42)) 

comes from the series expansions of the log function. We will show that another expressions which 
are less then (c) are positive in two steps. 

1. c > 1 : It is known from the calculus that if c > 1 then 

c- 1 



logc > 2- 



c+1 

Substituting 2^1 into p^(c) instead of logc one arrives at the expression 
1 

c 4 (l + c) 
which is positive if c > 1. 



(401c 5 + 185c 3 (c 2 - 1) + 93c 4 + 8c(c - 1) + 78) 



2. < c < 1 : The functions /i(c) = 8c 4 logc + 1, / 2 (c) = -(3c 3 + 2c 2 - c + 2) logc and 
/ 3 (c) = (444c 4 - 153c 3 - 50c 2 + 4c + 12) are positive on [0, 1]. From this follows that 

p W = 1 . (72/i(c) + 72/ 2 (c) + 8/ 3 (c) + 312) > if < c < 1. 



We proved that the function ti(c) (was defined by (|23|l 'l is concave. In the concavity proof of r 2 (c) 
(was defined by l|23Jl l we work with the function 

r*(c) =c 3 (c- l) 4 log 4 c-^(c), 

that is 

r*(c) =2c 3 (c + 2)log 4 c- 2(c- l)(3c 2 - 3c + l)log 3 c+ (3c- 2)(c- 3)(c- l) 2 log 2 c 
- (3c 2 - 12c + ll)(c - l) 3 logc - 9(c - l) 5 . 

To prove the concavity of t 2 (c) enough to show that t* (c) is positive. Since one can check the following 
limits 

lim r*(c) = lim r* (1) (c) = • •• = lim r* (5) (c) = 

c— >1 c — >1 c — >1 
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enough to show that the function 

p(c) = ^r* (5) (c) =48(2c- l)log 3 c+ -|(100c 4 - 11c 3 + 12c 2 + 12c + 12) log 2 c 

- ^|(90c 5 - 246c 4 - 216c 3 - 62c 2 - 9c + 78) logc 

- 4r(951c 5 - 642c 4 - 403c 3 - 36c 2 + 100c - 42) 
c J 

is positive if < c < 1 and negative if 1 < c. Using the previous methods again one can check that 

lim p(c) = lim p'(c) = lim p"(c) = 0. 

c— *1 c— *1 C— *1 

This means that enough to show that the function 

c 5 

— ■ p {2) (c) = (144c 4 + 72c 3 + 72c 2 + 216c + 432) log 2 c + (-180c 5 + 888c 4 - 78c 3 - 92c 2 
- 306c - 1440) logc - (1221c 5 - 1056c 4 + 282c 3 + 150c 2 + 273c - 870) 
is positive if < c < 1 and negative if 1 < c. Since 

limy ■pW(c)=0 
c— >1 z 

we will show that the function 

V[C) = 576c3 + 216c 2 + 144c + 216 ' ~dc [~2 ' 9 (c) ) (24) 

that is 

l2 2(450c 5 - 1920c 4 + 45c 3 + 20c 2 - 63c - 432) , 

^)=- l0 S c + c(576c3 + 216c 2 + 144c + 216) bgC 

6285c 5 - 5112c 4 + 924c 3 + 392c 2 + 579c + 1440 
+ c(576c 3 + 216c 2 + 144c + 216) 

strictly positive. 

We will show this inequality in two steps using an appropriate approximation of the log function. 

1. c > 1 : It is known that if c > 1 then y/c > logc. If one substitutes — c into the previous formula 
instead of — log 2 c then one decreases the 77(c) function. We decrease again the function and 
multiplying by 1152c(llc 3 + 2c + 3) we show that 

-(91905c 4 - 9081c 3 - 9064c 2 + 2016c + 13824) logc + (33c 3 + 6c + 9) (1855c 3 - 1776c 2 - 72) 

is positive. Since the coefficient of the log c is strictly negative if c > 1 enough to show that the 
function 

1 3(llc 3 + 2c + 3)(1855c 3 - 1776c 2 -72) 

(c) - - og c + c( - 91905c4 - 9081c3 - 9064c2 + 2016c + 13824) 

is positive if c > 1. One can check that h(l) = > and 

864(87213719c 4 + 366253c 3 + 339642c 2 + 218160c + 10368) + c 4 q(c) 



h'(c) 



c 2 (91905c 4 - 9081c 3 9064c 2 + 2016c + 13824) 2 
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where g(c) is a polynom of c. To see that h'(c) is positive if 1 < c enough to use the inequalities 

q(l) , , q^(l) , q^(l)>0 

and the inequality 

g (4) (c) = 2376(852418875c 2 - 482743225c - 4909728) > if 1 < c. 

2. < c < 1 : Since 77 ( 1 ) > (77 was defined by J23J) enough to show that 

36c 2 (8c 3 + 3c 2 + 2c + 3)-^<0. 

dc 

The previous expression is 

(3600c 8 - 1908c 7 - 6876c 6 - 5552c 5 - 20058c 4 + 12888c 3 + 3210c 2 + 1080c + 1296) • log c 
+ (28740c 8 + 4845c 7 + 2991c 6 + 22155c 5 - 35730c 4 - 25475c 3 - 7833c 2 - 3933c - 3456). 

After increasing this expression and dividing by c 3 we will show that the function 

-2c(954c 3 + 3438c 2 + 2776c + 10029) • log c + (28740c 5 + 4845c 4 + 2991c 3 

+ 22155c 2 - 35730c - 25475) 

is still negative if < c < 1. Since the coefficient of logc is negative if < c < 1 we will show 
that the function 

28740c 5 + 4845c 4 + 2991c 3 + 22155c 2 - 35730c - 25475 
Tj (c) - ogc 2c(954c 3 + 3438c 2 + 2776c + 10029) 

is positive if < c < 1. Since 77* (1) > enough to show that 

-2c 2 (954c 3 + 3438c 2 + 2776c + 10029) 2 • 

dc 

is positive. Computing the previous expression we get 

(472766109c 2 - 59724482c + 25548875) + a 3 c 3 + a 4 c 4 + • • • + a 8 c 8 
where 03, 0,4, . . . , as > and the first part is positive if < c < 1. 

□ 

The monotonicity of the function /32,fc(a, b) means that for given positive numbers 1 > a > b and 
eigenvalue 1 > Xk > 0, the function /3 2 ,fe (a — y,b + y) is strictly increasing in the variable x € [0, 3l ^\ ■ 

Using the a = j-, b = j- and c = notations the monotonicity means that the function 

/?2,fc(c + x, c — y) decreasing in the variable x € [0, c]. 
Let us define the following functions: 

Mu) = ^t— + t^— Mu) = ^— + ^—- (25) 
u log u 1 — u log u 1 — 11 

Using the equalities (|llf> 

Afc • (3 2 ,k(a,b) = 3^,3(2;, c) + 2qp{x 1 c) + 2rp(x,c), 
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where 



771 (c — X, C + X, l) 2 771 (c — X, 1, l)m(c + X. 1, 1) 

wa(x,c) = — ? ■ — — ; -r— -7— — , rp(x,c) = — — - — ■ -— , 

m(c — x, c + x)m(c — x, l)m(c + x, 1) m(c — x, l)m(c + x, 1) 

m(c — x, c — x, c + x)m(c — x, c — x, 1) m(c + x, c + x, c — x)m(c + x, c + x, 1) 



qp(x,c) = 



m(c — x, c — x)m{c — x, l)m(c — x, c + x) m(c + x, c + x)m(c + x, l)m(c — x, c + x) 
These function can be computed explicitly: 

log(c + x) c — x — 1 



tp(x,c) +tp(-x,c) , log(c-x) c + x-1 

wp(x,c) = — , where tp(x,c)- 



- 1 



x(log(x + c) — log(x — c)) ' 
0i (c + x) — (f>i(c — x) 



qp(x,c)= 9l ,p(x,c)+q 2lP (i:,c) where \ " logfaH- c) - log(aj - c) 

rp(x, c) = -0 2 (c - x)(j) 2 {c + x). 
The monotonicity of the function /3 2 ,k comes from the inequality 

— (3wp(x, c) + 2qp(x, c) + 2rp(x, c)) < 0. 

This inequality is rather difficult to prove but it can be separated into four seemingly monotone part. 
Theorem 4.3. If for all positive parameters c the following conditions hold for all x £E]0, c[ 

1: ^-{wp{x,c) +qp(x,c)) < 2: (2wp(x, c) + q h p(x, c) + rp(x, c)) < (26) 

then the (3 2 ,k part °f the scalar curvature is monotone, that is for given 1 > a > b positive numbers and 
1 > Afc > eigenvalue the function Pi y k(o, — x,b + x) is strictly increasing in the variable x £ [0, ^i^l ■ 



Lemma 4.1. The condition 4 fulfils in the previous theorem. 
Proof. The condition 

dx ' 



d -(rp(x,c)) < 



is equivalent to the condition 

<l>2( c + x ) 4>' 2 {c~x) 

4> 2 {c + x) <j) 2 (c-x)' 
where the function 4> 2 was defined by l|25|l . so enough to prove that 

&(«)Y 



<f>2 (u) 



> 



or equivalently 

p(u) = <h(u)<fc(u) ~ {^ 2 {u)f > 0. 
We will show that p*(u) = (1 — u) 4 log 4 zi • p(u) is a positive function. Since 

p4i) = pi 1 \i) = P i 2) (u) = pi 3 \i) = o 
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enough to show that the function t(m) = m 4 • pi 3 '' (u) that is 

t(u) = 2(-3u 3 - 6m 2 + w-6) log 2 u - 2(1 - u)(2u 3 - 9m 2 - 3m + 2) logu + 2(1 - u) 2 (6m 2 + 5m + 11) 
is negative if < u < 1 and positive if 1 < u. Since 

r(l) = r«(l) = p( 2 >(l) = t^(1) = r< 4 )(l) = 
enough to show that the function 

that is 

n(u) = (96m 4 - 72m 3 + 48u 2 + 8m + 144) log u - 8(1 - m)(61m 3 + 28m 2 + 28u + 30) 
is increasing or equivalently 

o 

tP(u) = (384u 3 - 216m 2 + 96u + 8) logw-(256u 4 - 108m 3 + 6u 2 + 3m + 18) > 0. 

M 

After dividing the previous inequality with the coefficient of logw which is strictly positive if u is 
positive, the inequality is the following: 

„ N , 16 18 1 2484m 2 - 1284m + 655 „ 

rf(u) = 1 ° SM+ y + - + 3- 48M 3 -27M 2 + 12M + l >0 - (27) 
We will show this inequality in three steps. 

1. < u < | : Since d (|) > enough to check that d*(u) = u 2 (48m 3 - 27m 2 + 12m + l) 2 • d'{u) 
is negative if < u < \. Since d*(0) < enough to show that d'^iu) < and this inequality 
follows from the equation 

d'„(u) = -m 5 (25920 - 16128m) - u 2 (-37245m 2 + 5400m 3 + 5000) - (5260m 2 - 2904u + 431), 

where d'^(u) is a sum of three negative terms. 

2. i < u < 1 : One can check that on this interval the function 

. . 1 2484m 2 - 1284u + 655 

v ; 3 48m 3 - 27m 2 + I2u + 1 

is concave hence 

c(u) > c(l)-c(V2) {u _ i/2) + c(i/2) if 1 < u < i. 

Substituting the right hand side of the previous inequality into the inequality Ij27(l we get that 

, , . 113323 18 13283 
10S(U) " - 2550- + u + ^2T U > °- 

Let decrease the left hand side of the previous inequality and show that 

18 1 

cAu) — log m — 45 H h 30m > if - < u < 1. 

M 2 

This function has only one stationary point which is a local minimum on this interval at mo = 
^M^i and Cjt ( Uf) ) > o. 
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3. 1 < u : Since d(l) > (was defined by O) enough to show that c(u) = w 2 (48u 3 - 27m 2 + 
12u + l) 2 • d'(u) is positive if 1 < u. Since 

c(l), C '(l), C ( 2 )(l) ;C ( 3 )(l), C W(l)>0 

enough to show that (u) > if 1 < U. This inequality follows from the equality 

c^)(u) = 5806080w 2 - 3110400u + 893880. 

□ 

The conditions 1,2,3 in Theorem l4.3l hopefullv can be proved similar way. The numerical evidences 
for previous theorem is the following: 

There was chosen 10 parameter c uniformly from the interval [0, 10 3 ] and for each chosen c the 
conditions 1,2,3 were numerically checked at 10 4 point uniformly from the interval [0,c]. This test 
can be view as a numerical evidence that the conditions 1,2,3 are valid. 

The summand j^l (a, b) in expression l|21|) can be computed using equalities (|11|) . I|13|) 

3w(a, X k ) + 3w(b, X k ) + q(d, X k ) + q(b, X k ) + r(a, X k ) + r(a, A fe ) + d(a, X k ) + d(a, X k ) 



7ki(a,b) 



A, 



where a = b = X k — 4^ and 

A; ' Al ' K A, 



w ( x , c ) = 7 — rn — ( 1 — ^ ~ — ~ ) ( + t — ~ ) > q( x > c ) = fai )^ (-) > 

(c — 1 ) log c \ log x — log c x — c ) \ log x 1 — X J \ c J 

r(x,c) = —02 (-) <f>2 ( - J , d(x,c) = <j>2{c)4>2{x). 

X \XJ \x J 

The function cf>i (c) and 4>2 (c) were defined by equation 1251) . The monotonicity of summand 7^ (o, 6) 
means that the function ^) k i{a — x,b + x) is increasing in the variable x € [0, ■ This property 
follows from the inequality 

d 2 

\3w{x, c) + q(x, c) + d(x, c) + r(x, c)J < 0. 

This inequality is rather difficult to prove but it seems that the left hand side of the previous inequality 
is the sum of two negative functions. 

Theorem 4.4. If for all c positive parameters the following conditions holds for all positive x 

d 2 d 2 
1: -j^(2w(x, c) +d(x, c)) < 2: (w(x, c) + q(x, c) + r(x, c)) < 

then the jk,l part of the scalar curvature is monotone. 

The conditions (1) and (2) in the previous theorem were tested numerically in the parameter space 
x, c G [0, 10 7 ], and the computation confirms the validity of the conditions. 

Note that equation (|21|l not the only reasonable decomposition of the scalar curvature formula 
(|15[) . There was tested many decomposition and the formula (|21|l seemed to be the most appropriate 
one. One can prove even more about summands in equation (|21|l : 
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1. a(a,b) (defined by fT7fr): 

For given parameters a > b > the functions 

if (a — x, a — x, b + x) + tp(b + x, b + x, a — x), 
if (a — x, b + x, a — x) + if(b + x, a — x, b + x) 

are strictly monotonously increasing if a; € [0, ■ From these follows Theorem (14.1(1 . 

2. Pi, k {a,b) (defined by fijfy): 

For given parameters 1 > a > b > and 1 > > the functions 

if (a — x, a — x, Afe) + ip(b + x, b + x, X k ), 

(f(a — x, a — x, Afe) + f(b + x, b + x, Afe) + if (a — x, Afe, a — x) + f{b + x, Afe, b + x) 
are strictly monotonously increasing but the function 

f(a — x, Afe, a — x) + f(b + x, X k} b + x) 
is not increasing if x S [0, s ^-\ ■ 

3. f3 2 , k {a,b) (defined by JTHJ)): 

For given parameters 1 > a > b > and 1 > Afe > the functions 

if (a — x, b + x, Afe) + f{b + x, a — x, Afe), 

p (a — x, b + x, Afe) + f(b + x, a — x, Afe) + ifi(a — x, Afe, b + x) + f(b + x, Afe, a — x) 
seem to be (numerically tested) strictly monotonously increasing but the function 

f(a — x, Afe, b + x) + f(b + x, Afe, a — x) 
is not increasing if x £ [0, ^i^l . 

4. 7fc,i(a,6) (defined by (HQ: 

For given parameters 1 > a > b > and 1 > Afe , A; > the function 

f{Xk, a ~ x, A/) + f{X kl b + x, Xi) 

seems to be (numerically tested) strictly monotonously increasing but the function 

if{a - x, Afe, A ; ) + f{b + x, Afe, A/) 

is not increasing if a; € [0, . 

If one defines another symmetric function 

7fe ; (a, b) =f{a, X k , Xi) + if(b, X k , A;) + f(X k ,a, A;) 

+ if(X k ,b, A/) + if(Xi, A fc , a) + if(Xi,X k ,b) 

then the function 7^ , (a — x, b + x) will be not increasing if the ratio of Afe and A; is large enough 
(- 15000). 
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5 Scalar curvature on the real and complex state spaces 

The real density matrices form a submanifold of the complex density matrices. The curvature tensors 
of a general submanifold can be very different from the curvature tensors of the manifold. For example 
one can think a circle (strictly positive scalar curvature) as submanifold of the plane (scalar curvature 
is 0). From this point of view one can reformulate Petz's conjecture for real density matrices. The 
question arises naturally: is there any connection between Petz's conjecture for real and complex 
matrices. 

Theorem 5.1. If the scalar curvature is monotone with respect to the majorisation on the space of 
complex density matrices then it is monotone on the space of real ones. 

Proof. From equations Q15[l. (|16|) we get the equation 

n 

Scal K (D) = - Scalp) + - £ v(X k ,Xi), 

k,l=l 

where the function v(Xk, Xi) was defined by (|13[1 . The monotonicity follows from the next two state- 
ments. 

1. For every eigenvalues A& > A; the function v(Xk — x,Xi + x) + v(Xi + x, Xk — x) monotone 
increasing in the variable x € [0, h 7 ] ■ 

2. For every eigenvalues Xk > A; and Xj the function v(Xk — x, Xj) + v(Xi + x, Xj) + v(Xj, Xk — x) + 
v(Xj, Xi + x) monotone increasing in the variable x € [0, ^ ] ■ 

Using equalities ^1 If) . I|13|l the function in the first statement can be written in the form 



v(X k - x, Xi + x) + v(Xi + x, X k - x) 



Xk + Xi 

where 



. . A fe - x . , 3 2c + 1 c + 2 

C(X) = , K(C) = 5 ; — h — -T7T. (28) 

Xi+x' w 2clog 2 c c(c- 1) log c 2(c-l) 2 

Since the function c{x) is decreasing enough to show that the function 



d(c) = (1 + c)k{c) + ( 1 + ~ K (\ 

is decreasing if c > 1. This will follow from the negativity of the function 

4(c) = (c- l) 3 log 3 c-d'(c). 

Since 



enough to show that 



lim 4(c) = lim<(c) = 

C— >1 



T(C) = (l-c)(6^ + 18c 2 + 6) ' <(C) < ° if 1 < C 



Since lim c _+i r(c) = enough to show that the function 

3c(c 4 + 3c 2 + l) 2 _ t =(6c6 + 22c5 + 5gc4 + 3gc 3 + 59c2 + 22c+ 6)bgc 
(1 - c) 2 

+ ^(l-c 2 )(c 4 + 6c 3 + 6c+l) 
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is positive. The coefficient of the fog c is positive, so enough to show that 



5„ 2 . c 4 + 6c 3 + 6c+l . . 

(c) = logc+ -(1 - c 2 ) ■ <0 ifc>l. 

W 8 2 V ' 59c 2 + 22c + 6 



Since lim c _ + i r*(c) = enough to show that r^(c) is positive and this follows from equation 

, , . a + a\c + a 2 c 2 H h ai 2 c 12 



c(6c 6 + 22c 5 + 59c 4 + 36c 3 + 59c 2 + 22c + 6) ' 

where ao, . . . , d\2 are strictly positive numbers. 
This completes the proof of the first statement. 

Using the equalities 1)11(1 , 1(1 3(1 the function in the second statement can be written in the following 
form 

v(X k - x, Xj) + v(\i + x, Xj) + v(Xj,X k —x) + v(Xj,Xi + x) = -^-(n(Xk - x) + n(Xi + x) 

+ p(X k -x) + p[Xi + £), 
where Xk — ^,Xi = y-,x = j-, the k(c) function was defined by ((28(1 and 

3 c + 2 1 + 2c 



Pic) 



21og 2 c (c-l)logc 2(1 -c) 2 



Using the same arguments as in the proof of Theorem l4.2l enough to show that the function k(c) + p(c) 
is concave for all positive c. Let us define the function 

d(c)= {c ~^ g4c ■( K "(c)+p"(c)). 

To prove the second statement enough to show that d(c) < for every positive c. Since 
lim d(c) = lim d'(c) = lim d (2) (c) = ■ • • = lim d (6) (c) = 

c — >1 c — >1 c — >1 c — >1 

enough to show that 

24 

r(c) = d (6) (c) =1440 log 3 c - ^(30c 6 - 211c 5 - 198c 4 + 207c 3 + 18c 2 + 54c + 90) log 2 c 

& 

+ 4(1 - c)(351c 5 + 1306c 4 - 29c 3 + 1120c 2 + 957c + 765) logc (29) 
c° 

+ 4(1 - c) 2 (1249c 4 + 1132c 3 - 744c 2 - 872c- 705) < if < c. 

c b 

Note that 

lim t(c) = lim t'(c) = 0, lim(cV(c))' < and r(4), r'(4), (cV(c))'(4) < 0. 

c — >1 c — >1 c — >1 

To prove the inequality (|29|1 for < c < 1 and 4 < c parameters enough to check that the function 

Mc) = (c 6 -T>(c)y (so) 

is positive if < c < 1 and negative if 4 < c. We will prove the inequality 1(29(1 in four steps. 
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1. < c < | : After decreasing the function t»(c) and dividing by c 2 one arrives at the inequality 

2 _4483c 2 + 108, _ 76944 16208 56520 
?/>(c) = -78624 log 2 c + 96 5 log c - 53008c + 194856 — + — — . 

C A c c z c 

Since ijj (j) =0 enough to show that 

cV(c) 



V*(c) 



157248c 3 + 430368c 2 + 31104 
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logc+ —(3313c 5 - 31707c 3 - 2026c 2 - 648c + 14130) > 0. 
c 



Since V 7 * (5) > enough to show that 

p(c) = 3c 2 (1638c 3 + 4483c 2 + 324) 2 ^(c) < 0. 
The function p(c) is a polynom and can be checked that p(c) <0if0<c<^. 

2. i < c < 1 : After decreasing the function r*(c) (which was defined by one arrives at the 

inequality 

-819c 2 log 2 c + 4483 ° 2 + 108 log c - 553c 3 + 2029c 2 - 800c - 169 + ^ > if - < c < 1 
c c^ 2 

From the Taylor-expansion of the functions one can check that 
4483c 2 + 108 



c 

where 



■ log > a + 01 (c - 1/2) + a 2 (c - 1/2) 2 , -819c 2 log^ c> fe + 61 (c - 1/2), 



4915 W2 

a = —2-, ai = 4915 -4051 log 2, a 2 = 3187 - 864 log 2, 

bo = _ 8191og 2 2 ^ 6 1 = 8191og2(l-log2). 

After substituting the logc functions with the polynom of (c— 1/2) we get the following inequality 

^ (2212c 5 + a 4 c 4 + a 3 c 3 + a 2 c 2 - 2352) > 
4c^ 

where a 4 = 3456 log 2 - 20864, a 3 = 3276 log 2 2 + 9472 log 2 - 3712 and a 2 = -819 log 2 2 + 
4230 log 2 + 7319. One can check that the polynom in the parenthesis is strictly negative if 
\ <c< 1. 

3. 1 < c < 4 : Let increase the function r(c) (defined by (|23fl ): substitute — 2500(c— 1) instead of 
24 

-(30c 6 - 211c 5 - 198c 4 + 207c 3 + 18c 2 + 54c + 90) log 2 c, 



and -2000(c- l) 2 instead of 
8 



-r(l - c)(351c 5 + 1306c 4 - 29c 3 + 1120c 2 + 957c + 765) logc 
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and 4r • (1 - c) 2 (1249c 4 + 1132c 3 ) instead of 

^5 (1 - c) 2 (1249c 4 + 1132c 3 - 744c 2 - 872c - 705). 

That these substitutions increases the function r(c) can be check using the Taylor-expansion. 
The new inequality is 

ri (c) = 1440 log 3 c+ 4 ^7 ^ • (500c 3 - 1124c 2 + 117c + 1132) < if 1< c < 4. 
Let increase the function n (c) : 

73, (c) = ri (c) + 2000 (C ~ 1)( 2 C ~ 2)2 . 

c z 

Since lim c ^i 72(c) = to show that 72(c) is negative enough to check that ct' 2 (c) < 0. Let 
increase the function ct' 2 (c) < 

73 = ct^c) + 2000 (c- + 8(c-4) 2 -72. 

Since lim c _>i 73(c) — to show that 73(c) is negative enough to check that 

> f y 4(-996c 4 + 608c 3 + 2985c - 3472) n 
rg(c) = 8640 logc+ — — -<0 if 1 < c < 4. (31) 

If 2 < c < 4 then one can substitute c — 1 instead of log c and one can check that the inequality 

-996c 4 + 2768c 3 - 2160c 2 + 2985c - 3472 .„ „ 
- if 2 < c < 4 

c z 

holds. If 1 < c < 2 then one can check the inequality H31J) holds using the inequality 
c-1 (c-1) 2 (c-1) 3 (c-1) 4 (c-1) 5 (c-1) 6 



2c 2 3c 3 4c 4 5c 5 



> logc if 1 < c < 2. 



4. 4 < c : After increasing the function r*(c) (which was defined by H30(| s | one arrives at the 
inequality 

(86400c 3 + 59616c + 2592) log 2 c + —(-2985c 5 + 5840c 4 + 3126c 2 + 216) logc 

c 

+ 4( 21 84c 6 + 24357c 4 + 204c + 7065) < if 4 < c. 

The coefficients of log 2 c and log c are positive and c — 2 > log 2 c, c — 2 > log c if 4 < c therefore 
one way to increase the function is to substitute c — 2 instead of log 2 c and log c. Then the new 
inequality is 

c 4 (5970c 3 - 27948c 2 + 30560c - 16855) + (17364c 3 - 216c 2 + 796c - 2355) > if 4 < c. 

This inequality is holds since the left hand side is the sum of two positive functions on the 4 < c 
interval. 

□ 
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6 Conclusions 



We showed that Petz's monotonicity conjecture for the scalar curvature of the Kubo-Mori metric 
follows from more elementary inequalities. The key idea was to find a good grouping of summands in 
the expression of the scalar curvature. The proof of some inequality was given using elementary but 
brutal computations. We proved that if Petz's conjecture holds for the manifold of complex density 
matrices then it is also true for the manifold of real ones. 

The higher order derivatives of the summands of the scalar curvature play central role in these 
computations. The scalar curvature is a complicated expression of the Kubo-Mori metric, which can 
be derived from von Neumann entropy. It seems that these higher order derivatives (of the entropy 
function) responsible for the monotonicity of the scalar curvature. It would be good to find a proof 
for Petz's conjecture which is based only the easily computable properties of entropy function. 

Acknowledgement. I would like to thank Dr. Denes Petz for stimulating discussions and 
valuable remarks. This work was partially supported by OTKA32374. 
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